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ABSTRACT 

I revisit the Ceph eid-distance determ i natio n to the nearby spiral galaxy Ml01 (Pin- 
wheel Galaxy) of IShappee fo Stanekl (120111 ). in light of several recent investigations 
questioning the shape of the interstellar extinction curve at A « 8,000 A (i.e. I- 
band). I find that t he relatively steep extinction ratio Aj/E(V — I) = 1.1450 
( Fitzpatrick fe Massal l2007h is slightly favoured relative to Aj / E(V — I) = 1.2899 
( Fitzpatriddll999l ) and sig nificantly favoured re lative the historically canonical value 
of Aj/E(V — I) = 1.4695 (ICardelli et al.l[l989l) . The steeper extinction curves, with 
lower values of Ai/E(V — /), yield fits with reduced scatter, metallicity-dependences 
to the dereddened Cepheid luminosities that are closer to values inferred in the lo¬ 
cal group, and that are less sensitive to the choice of reddening cut imposed in the 
sample selection. The increase in distance modulus to M101 when using the preferred 
extinction curve is A/i ~ 0.06 mag, resulting in an estimate of the distance modulus 
to M101 relative to the LMC of A/iLMC ~ 10.72 ±0.03 (stat). The best-fit metallicity- 
dependence is dMj/d[ O/H] ss (—0.38± 0.14 (stat)) mag dex -1 . 

Key words: distance scale - stars: variables: Cepheids - ISM: dust, extinction 
galaxies: individual: (M101) 


1 INTRODUCTION 

Interstellar extinction causes luminosity sources to appear 
fainter and redder, and is thus somewhat degenerate with 
distance which causes luminosity sources to appear fainter 
but not redder. In principle, this difference (the colour offset) 
can allow one to break the degeneracy between interstellar 
extinction and distance. If one assumes a consistent total-to- 
selective extinction ratio Ri = Ai/E(V — I) = Ai /(Ay — 
Ai), one can estimate a “WesenheilQ magnitude” (iMadord 
1 1982h that is nominally reddening-independent: 

Wi = I — Rix(V — I). (1) 


The value s of Ri regularly us ed in the literature range 
from 1.60 i Paczvnski et al.lll994l). 1.55 rtMaiaess et ~al] 201ll : 
Ngeowll2012l). down to 1 . 45 (iMacri et all 20061 : Gerke et all 


201 ll : Shappee fe Staneklfeoill) . There is certainly a stron g 


empirical basis for this assumption. ICardelli et al.l (1 1989h . 
in their seminal investigation of extinction from the ultra¬ 
violet to the near-infrared, yielded parametric fits for the 
extinction curve as a function of wavelength and the free- 
parameter Rv = Ay/E(B — V). For the case Ry = 3.1, 
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their f unctions yield Ri ~ 1.45. The results of lCardelli et al.l 
(I 19891) are a linchpin of astronomy as a whole, with over 
5,000 refereed citations at the time of writing. Thes e values 
have h ad convincing independent support elsewhere. [Udalskil 
(l2003l) measured Ri = 1.44 ± 0.03 toward the Large Magel¬ 
lanic Cloud (LM C), which is the ancho r of the extragalac- 
tic distance scale (|Freedni;ML_eLjdJ|200l|) a nd m any investi¬ 
gations of Cepheids. IPeicha fe Kochanekl ll2012l ) developed 
a comprehensive physical model and matched to ~177,000 
photometric measurements from a sample of 287 LMC, 
Small Magellanic Cloud (SMC), and Milky Way Cepheids. 
They measured a best-fit mean exti nction curve o f Ry ~ 
3.127, ass uming the forma l ism f rom ICardelli et all (Il989l) . 
Similarly. IPeicha fe Prietol (120141 1 . matched to ~6,800 pho¬ 
tometric measurements over a broad range of wavelengths, 
and showed that t he Ry = 3.1 extinction curve from 
ICardelli et al.l dl989l) is an effective fit to the UBVRI light 
curves of Type-II plateau supernovae. 

There has been recent evidence that the extinction 
curve of even the diffuse interstellar medium may be vari¬ 
able on a syst ematic basis, or sim pl y distinct than as pa- 
ramet erised by lCardelli et all ll 198911 . ISchlaflv fe Finkbeinerl 
d201ll) studied reddening values toward Milk y Way halo 
stars, and found that the exitnction curve of iFitzpatrickl 
d 19991) . with its different method of derivation with the ex- 
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tinction curve fit to some cubic spline anchor points rather 
than seventh degree polynomials among other distinctions, 
and its steeper extinction in the wavelength range 6,000 
A sj. A $^9,000 A, was mo re consistent with observations that 
of ICardelli et al.1 (Il989lf . This finding contributed to a re - 
vision of the SFD extinction m aps (ISchleeel et al.l 1 1998h . 
Similarly, ISung fe Besselll (l2014l f showed that reddening to¬ 
ward open c lusters with i n 3 Kpc of the Sun is well-fit by 
an Rv = 3.1 [Fitzpatrickl (il999l l extinction curve, which pre¬ 
dicts Ri ~ 1.27 for hot stars. Further, it is by now well- 
documented that the inner-parts of the Milky Way, which 
are rigorously investigated, are fit by a steepe r standard ex¬ 
t inctio n curve than the Rv = 3.1 curve from lCardelli et al.l 
(1 198911 . In their study of 17,000 Galactic bulge RR Lyrae 
stars i dentified in OGLE-I I I pho tometry (ISzvmanski et all 
l201lli . Pietrukowicz et all (120121 ') measure Ri = 1.080 ± 
0.007. iNataf et al.l ( 2013l l used red clump stars (which are 
~2,000 Kelvin colder and ~ 10x more metal-rich) mea¬ 
sured with the same OGLE-III photometry and obtained 
< Ri >= 1.218, with significant variations confirmed span¬ 
ning a range 1.0 < Ri < 1.4. The large shifts in the extinc¬ 
tion curve are observed to take place toward sightlines closer 
to one another than 0.5 degrees on the sky, or a transverse 
distance no g reater than 75 parsecs. Finally, from the extinc¬ 
tion study of I Berdnikov et ahl (1 19961 '). which based purely off 
photometry of Cepheidfl, Ri « 1.07. 

Reading this sequence of findings some may argue that 
the astronomical com munity should simply switch to using 
thejFitzjoatrick ( 199j|) extinction curve as default over that 
of lCardelli et ali ( 1989i l. That may turn out to be the correct 
course of action, but as of now there at least three causes for 
concern with this approach. First, as discussed earlie r in the 
introduction, the extinction curve of lCardelli et al.l (1 19891) 
does have convincing independent empirical support, and 
this wo uld need to be explai ned if the Rv =3.1 extinction 
curve of lCardelli et al.l l|l989l f turns out to be incorrect. Sec¬ 
ond, the evid ence is curre ntly overwhelming that the extinc- 
tion curve o f |Fi tzpatricl3 d 199914 fails in the near-IR, where 
iFitzpatrid^ ill 99(1) predicts Aj/A k —— 2.31 that is nearly 
independent of Rv- However, the most precisely and accu¬ 
rately measured value in the Galactic as tronomy literature is 
Aj/Ak ~ 3.02 llNishivama et al.l 120091 4 . Finally, buoyed by 
more abundant data and a different mode of analysis, Fitz¬ 
patrick hims elf argued for a differ ent mean Galactic extinc¬ 
tion curve in lFitzpatrick fc Massal J2007 ) (se e th eir Table 5), 
which is as di stinct from IFit zpatric kl J 19991 ') as iFitzpatrickl 
(Il999l4 is from ICardelli et al.l ( 1989l 4~^ 

Separately and independently to these issues, evidence 
has recently emerged that the wavelength-dependence of in¬ 
terstellar extinction is a function of the s pectral type and 


the sp ecific-star-formation rate of a galaxy (Kriek_&_Conro^ 


2013), and even location within a galaxy (IZasowski et ah 
200£ ; (Phillips et al.ll2013l : [Dong_et_al. 201 A L Durin g the ref¬ 


eree process for this paper, Fausnaueh et al.l (120141 4 found a 
best-fit extinction curve toward the Cepheids of NGC 4258 


2 From Table 4 of I Berdnikov et ahl jl996l i one can derive Rj = 
0.84, but that is because the /-band filter used in that study is 
centred at A Ri 8, 800A, which is at a longer wavelength than the 
standard Landolt filter. 


of Rv = 4.9, which is extremely shallow, using a combina¬ 
tion of UBVR photometry. 

Having worked extensively on some of these issues, 
I cannot argue at this time that there is a convincing 
framework unifying these disparate, seemingly mutually- 
inconsistent findings. What I am confident in, is that more 
investigation is needed, both of the shape of the extinction 
curve, how it varies with direction, and how this hitherto 
largely ignored source of systematic error has metastasised 
into the general astronomy literature and its library of ac¬ 
cepted values. There is an abundance of data available in 
a broad range of bandpasses for a broad range of objects 
observed in a broad range of contexts, and thus further re¬ 
search that could elucidate these issues is possible. 

Thus, I tackle the Cepheid distance to M101 for two 
reasons. First, as a probe of the interstel lar extinction curve 
in and of itself. The data products of IShappee fc Stanekl 
(l201lll are of high quality, with the photometry obtained 
from the Hubble Space Telescope (HST) and reduced via 
a suite of state-of-the-art software tools. The observations 
are of a spiral galaxy of similar mass and morphology as the 
Milky Way. As M101 is viewed nearly face-on, the reddening 
values are sufficiently small that a large numbe r of Cepheids 
can be identified. Indeed, the catalogue of lShappee fc Stanekl 
d201ll4 is one of the largest homogenously-measured extra- 
galactic catalogues of Cepheids, containing 1,227 Cepheid 
candidates. Second, this investigation can help discern what 
impact these uncertainties might have on the field of the 
extragalactic distance ladder , where Cepheids are a widely- 
used standardisable candle llMacri et al.l 120061 : iRiess et al.l 
l201ll ; iGerke et al.ll201 ll : lEfstathioull2014l l. As the era of pre¬ 
cision cosmology is now over a decade old dSpergel et al.l 
1200314 . with the reported error b ars growing smaller and po¬ 
tenti al tensions growing greater (|Planck Collaboration et al.l 
1201314 . the requirements for measurement precision are now 
higher-than-ever, necessitating diligent investigations of pa¬ 
rameters such as extinction coefficients. 


2 DEFINING THE PROBLEM AND 
METHODOLOGY 

As discussed in the introduction, virtually any extinction co¬ 
efficient in the range 1.05 < Ri < 1.45 can be assigned some 
empirical support ascribing it the designation of “the stan¬ 
dard extinction ratio for the diffuse interstellar medium.” 
That is problematic and a source of major systematic error 
previously undiagnosed and unquantified in the literature on 
the use of Cepheids to probe the extragalactic distance lad¬ 
der. In this paper, I test for three differ e nt cas e s. These are 
the ex tincti on curves of |Cardelli_et_alJ (1 198914 , IFitzpatrickl 
ljl999l4 , and iFitzpatrick fc Massal d2007l4 , which are shown 
in Figure [T] I c onvolve these c urves with the transmission 
functions of the lLandoltl (1 19921 1 V and I filters (also shown 
in Figure [TJ, and a synthetic spectrum of a T e fj = 5,750 
Kelvin, [Fe/H]= 0, logg = 1.0 star kindly sent to me by 
Luca Casagrande, where the pa rameters are chosen t o typ i- 
cal of Cepheids as modelled by lPeicha fc Kochanekl ll2012ll . 
The predicted extinction coefficients are listed in Table Q] 
At the request of the referee, I also list in Table [5] what 
the extinction coefficients would be for different assumed 
parameters. Shifting the temperature by 1000 Kelvin only 
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Table 1 . Predicted extinction coefficients for the Lan- 
dolt, 2MASS, ACS, and WFC3 filter transmission 
curves convolved with the synthetic spectrum of a T e g = 
5,750 Kelvin, [Fe/H]= 0, log g = 1.0 star as well as 0.40 

mag of extinction at A = 5500A. These are the coeffi¬ 
cients assumed by this study. 



CCM89 

F99 

FM07 

Aj/Eiv - 1) 

1.4695 

1.2899 

1.1450 

Au /A5500A 

1.5492 

1.5532 

1.5276 

Ab/A 5500 a 

1.2893 

1.3022 

1.2774 

Av /AI5500A 

1.0020 

1.0015 

0.9965 

Ar/A 5500A 

0.8187 

0.7854 

0.7683 

Ai / A5500A 

0.5963 

0.5641 

0.5319 

-Aj /^s^nnA 

0.2981 

0.2694 

0.2380 

Ah/A 55 00^ 

0.1825 

0.1716 

0.1410 

Aif s /A R5nn ^ 

0.1178 

0.1160 

0.0853 

Af5551VACs/A 5500 a 

1.0367 

1.0423 

1.0356 

^-F606VT,ACsM 5 500A 

0.9339 

0.9178 

0.9058 

Af 814(T, ACS / Assnn A 

0.5699 

0.5406 

0.5079 

Af.555 W, WFC3 / A 5500 ^ 

1.0482 

1.0521 

1.0414 

Af606VT,WFC3 / ^5500A 

0.9304 

0.9138 

0.9019 

AF814W,WFC3M<^nn A 

0.6001 

0.5671 

0.5354 

^4-F110VT,WFC3M 5 500A 

0.3379 

0.3184 

0.2855 

AF125VyWFC3/A 5500 y 

0.2881 

0.2686 

0.2372 

Af 160VC WFC3 / A RRnn A 

0.2042 

0.1907 

0.1603 


shifts the extinction coefficients by ~1% or less. The dom¬ 
inant uncertainty lies in the assumed extinction curve, and 
not the convolution of the stellar spectrum with which the 
transmission filters are convolved. Nevertheless, should one 
actually be attempting “1% cosmology”, then one should 
compute slightly different extinction coefficients for every 
Cepheid. With the two bandpasses used in this work, the 
shift in distance modulus is only ~0.003 mag, with the val¬ 
ues of Table 1 yielding the marginally smaller distances. 

I then follo w the sample select i on ou tlined in Section 4 
and Section 5 of Shappee fc Stanekl (1201 lit , which is based off 
the work of iMacri et al.l ( 2006|), including the same period- 
luminosity relations from Udalski et al. i a and sample- 
cuts based on colour and variability amplitude ratios meant 
to reduce blends. However, 1 make a few modifications to 
the methodology: 


(i) I use recursive 3-<r outlier rejection rather than itera¬ 
tive weighted-median sigma clipping to remove outliers. This 
consistently leads to 3 of the 231 clean Cepheids being re¬ 
moved from the fit, independent of the assumptions made 
in the fit. 

(ii) The errors in the relations are computed assum¬ 
ing homoscedastic errors on the data poin ts. I arg ue that 
the y 2 estimates from the m ethod of iMacri et al.l (l2006h 
and lShappee fe Stanekl (|201lh are problematic. The various 
sources of systematic error (such as the uncertain extinc¬ 
tion curve) are in need of characterisation, and these will 
add to the existing errors in unexpected ways, rather than 
simply as a re-scaling of existing errors. Further, there is no 
correlation given for the errors on various variables, such as 
the error on the mean 7-band magnitude or U-band mag¬ 
nitude of the Cepheids, which could be correlated. Thus it 
is actually impossible to compute a y 2 from the available 


Table 2. Predicted extinction coefficients for the Lan- 
dolt, 2MASS, ACS, and WFC3 filter transmission 
curves convolved with the synthetic spectrum of a T e g = 
6, 750 Kelvin, [Fe/H]= 0, log g = 1.0 star as well as 0.40 

mag of extinction at A = 5500A. 



CCM89 

F99 

FM07 

Aj/EiV-I) 

1.4524 

1.2721 

1.1304 

Au/A 55 ooa 

1.5395 

1.5407 

1.5151 

Ab/A 5500 ^ 

1.3129 

1.3233 

1.2979 

Av /^5500A 

1.0078 

1.0085 

1.0035 

Ar/A 5500 ^ 

0.8258 

0.7933 

0.7768 

Ai M5500A 

0.5969 

0.5646 

0.5325 

Aj A^5500A 

0.2900 

0.2702 

0.2388 

Ah M5500A 

0.1826 

0.1717 

0.1411 

Aks/A^ oo^ 

0.1178 

0.1160 

0.0853 

Af555VFACsM5500A 

1.0435 

1.0503 

1.0432 

^F606VF,ACSM 5 500A 

0.9491 

0.9360 

0.9244 

Af814W, ACS / A e Rnn a 

0.5711 

0.5417 

0.5090 

^F5551T,WFC3 / ^5500A 

1.0636 

1.0698 

1.0583 

^F6064V,WFC3 M5500A 

0.9447 

0.9309 

0.9195 

AF8141T,WFC3/A r ^^ m £ 

0.6016 

0.5684 

0.5368 

Afi 10 WVWFC3 / A 5500 a 

0.3429 

0.3234 

0.2904 

AF125TV,WFC3/A RRm A 

0.2895 

0.2699 

0.2385 

AF160WVWFC3/A RRnn A 

0.2045 

0.1910 

0.1606 


information even if one assumes that all sources of error are 
known, since y 2 is a correlation-dependent quantity. 

(iii) I fit for the metallicity-dependence of the dered- 
dened Cepheid luminosities concurrently with the fit for 
the relative distance modulus of M101 with respect to 
the L MC using a mult i linear least-squares fit. In con¬ 
trast, [Shappee fe Stanekl (j201ll l computed the metallicity- 
dependence using two sequential least-squares fit, with 
the distance modulus computed first independently of the 
met allicity-dependence, and the scatter to the first relation 
then r egressed with respect to th e metallicity. 

(iv) IShappee fc Stanekl (1201 lli set their metallicity zero- 
point to [0/H]= 8.50, which is thus their assumed metal¬ 
licity for the LMC. I instead use [0/H]= 8.36, as ar¬ 
gued b y lBresolinl (1201 if) . which was assumed bv lRiess et al.l 
(1201 ill and is consistent with the measurement of [0/H] = 
8.37 in HII regions and supernovae remnants of the LMC 
llRussell fe Dopitalll990l l. 

(v) As this investigation is concerned with quantifying 
one particular systematic error, the estimates of various pa¬ 
rameters state only the maximum -likelihoood values and the 
statistical errors. In contrast to IShappee fe Stanekl (l201lh 
who report estimates of the total systematic error separate 
from their statistical errors. 


3 RESULTS 

3.1 Metallicity-Dependence as a Free Parameter 

Assuming the extinction curve of ICardelli et al] (ll989i l. I 
obtain: 

A/mMC = (10.663 ± 0.031) + (-0.433 ± 0.145)([O/H] - 8.36) 

5 = 0.1680, 

( 2 ) 
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Figure l. The exti nction curves of ICardelli et al.l (I1989T) as the dotted-blue line, I Fitz Patrick! (Il999h as the short- 
dashed-green, andjFjtapa|ridO'£ Massal ( 200711 as the long-dashed-salmon line, as a function of wavelength. Also 
shown are the lLandoltl ( 1992(1 V and Landolt I filters in black. The different trends of the extinction curves over this 
wavelength regime lead to different predictions for the extinction ration Ai/E(V — I). 


where 5 is the 1-cr scatter on t he re lation. If I assume the 
extinction curve of lFitzpatrickl (Il999l l: 

A^lmc = (10.692 ± 0.031) + (-0.406 ± 0.144)([O/H] - 8.36) 

<5 = 0.1671, 

(3) 

and for the extinction curve of lFitzpatrick fe Massal (i2007lf : 
A^lmc = (10.716 ± 0.031) + (-0.384 ± 0.144)([O/H] - 8.36) 

S = 0.1674, 

(4) 

I reiterate that the errors quoted are random errors resulting 
from the assumption of homoscedastic errors on the data 
points. A complete inventory of the systematics involved in 
Cepheid determinations and their associated correlations is 
not available at this time and is beyond the scope of this 
work. Only one systematic is being focused on: the shape of 
the extinction curve. 

The relative distance estima te from Equations 2,3,4 is 
~0.03 mag greater than that from lShappee fe Stanekl ([2011) 
in spite of assuming the same extinction curve and nearly 
the same selection function. The offset is largely due to the 
downward adjustment of the assumed LMC metallicity by 
~0.14 dex. 

For the distance relations fit by Equation [5] through to 
Equation [ 6 ] in the following subsection, the same three out¬ 
liers and no others are always removed: F2-3598, Fl-3546, 
and F2-1458. Their respective inferred distances are ~0.60 
mag, ~-0.70 mag, and ~1.20 mag further than the relations. 
Due to these large offsets and the fact that it’s the same 
three outliers removed each time, I argue that the recursive 
3-cr outlier rejection does not detrimentally bias the results. 

The metallicity dependence derived, 7 = 

dMi/d[M/H] ~ 0.40 mag dex -1 , is substantially lower 
than the values of 0.72 mag dex -1 and 0.80 mag dex -1 


that IShappee fe Starred (|201ll i derive using two different 
methods. That may be due to the different method of 
fitting. The value derived here is closer to other literature 
estimates derived from other estimates, as will be shown in 
the subsequent subsection. 


3.2 Metallicity-Dependence Fixed to Literature 
Values 

It is of interest to verify how results differ when 
the metalli c ity -dependence is fixed to literature values. 

I Sakai et alj (|2004l f estimated 7 = —0.24 ± 0.05 mag dex -1 
by comparing tip of the red giant branch magnitudes to 
Cepheid magnitudes for 7 nearby galaxies images in V and 
I with HST , as well a s 10 a dditional values found in the 
literature. IStorm et al.l (1201 il l estimated 7 = —0.23 ± 0.10 
mag dex -1 for a sample of 51 Cepheids from the SMC, the 
LMC, and the Milky Way. I adopt the value of 7 = —0.23 
mag dex -1 though the two values are so similar that the 
effect of this choice is negligible. 

The relation o btained for the extinction curve of 
ICardelli et al.l Il989lj is: 

A/xlmc = (10.622 ± 0 . 011 ) 

S = 0.1687, ^ 

for the extinction curve of lFitzpatrickl [19991 1: 

A/xlmc = (10.657 ± 0 . 011 ) 

( 6 ) 

<5 = 0.1676, 

and for the extinction curve of lFitzpatrick fe Massal (l2007tl : 

A/tlmc = (10.685 ± 0 . 011 ) 

(5 = 0.1678, 

M101 is shifted to a distance that is ~0.04 mag further 
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Table 3. The slope of the shift in inferred distance mod¬ 
ulus per f magnitude increase in the reddening cut, as 
a function of assume d extinction curve. The curve of 
iFitzpatrick fc Massal 2007T) yields the smallest depen¬ 
dence to the choice of colour cut. 


CCM89 

-0.063 

F99 

-0.027 

FM07 

-0.017 


away in all three cases relative to the values obtained when 
the metallicity-dependence is floated as a free parameter. 
That offset comes from the product of the mean [O/H] of 
M101 relative to the LMC (~0.20 dex) and the shift in the 
slope of the metallicity-dependence (~0.15-0.20 mag dex -1 ). 

Though this was interesting as an exercise, I note that 
this assumption of a metallicity-dependence form the liter¬ 
ature is inconsistent with other assumptions in this inves¬ 
tigation. IStorm et al.l (120111) assumed a universal, constan t 
extinction function of Rv = 3.23 from lCardelli et aD (1 19891) . 
corresponding to Ai/E(V — I) tn 1.50. That’s marginally 
shallower than the LMC value of Ai/E(V — I) = 1.44±0.03 
dUdalskill2003l ). and substantially shallower than the Milky 
Way value o f Ai/E(V — I) ~ 1.27 suggested by the work of 
ISchlaflv fc Finkbeineij (1201 ll) . It may also be distinct from 
the SMC value. There is an urgent need for an investigation 
of the metallicity-dependence of Cepheid luminosities that 
accounts for extinction curve uncertainties and variations in 
the Milky Way Galaxy and its satellites. 
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Figure 2. TOP: The distance modulus of M101 rela¬ 
tive to that of the LMC as a function of the extinction 
curves and as the maximum reddening E(V — I )Max 
of Cepheids i ncluded in the calcu lation. The extinc¬ 
tion curve of iFitzpatrick fc Massal (120071 ) is the least 
se nsitive to the choice of E(V — I) Max, whereas that 
of ICardelli et al.l ( I1989T ) is the most sensitive. BOT¬ 
TOM: Same as top panel, but with the distance mod¬ 
ulus for E(V — I]m,v = 0.39 (the value assumed by 
IShappee fc Stanekll201ll) . shifted to zero. 


3.3 Shifting the Maximum-Colour Cut 

IShappee fc Stane~kl d 201 ll ) removed Cepheids with E(V — 
I) > 0.39 from their fit, on the basis that they might either 
be red blends, or have high-reddening values that would lead 
to an error linearly proportional to the error in the assumed 
extinction curve. Of interest would be to test which of the 
three extinction curves yields a distance least sensitive to 
the choice of E(V — I) cut (hereafter: E(V — /)mhx). For 
this calculation, I fix the metallicity-d ependence to the lit¬ 
erature value of 7 = —0.23 mag dex -1 dStorm et al.ll201ll) to 
restrict the number of degrees of freedom, though it is noted 
that the results are nearly identical if I fix the metallicity 
dependence to the values found in Equations [2] [3] and [4] 

With this test, of which I plot the res ults in Figure 
[2] the extinction curve of lFitzpatrick fc Massal (120071 ) is the 
m ost successful a t matc hing the data for M101, whereas that 
of lCardelli et al] (1 19891) which is the most standard in the 
Cepheid literature is the least successful. The inferred rel¬ 
ative distance modulus decreases rapidly as more reddened 
sources are in cluded when I assume the extinction curve of 
ICardelli et al.1 (Il989l ). exactly what one would expect if the 
extinction curve is too shallow. The distance mod ulus also 
decreases when I assume the extinc tion curves of lFitzpatrickl 
(119991 ) IFitzpatrick fc Massal (120071 ), but by a much smaller 
amount, ~ —0.027 and ~ —0.017 respective decreases in dis¬ 
tance modulus per 1 magnitude increase in E(V — Thyiax. In 
contrast, use of the extinction curve of lCardelli et al.l (119891) 
yields a much steeper sensitivity of ~ —0.063. These results 
are listed in Table [3]. 

The same conclusions are reached if one would rather 


Table 4. Scatter in distance modulus as a func¬ 
tion of colour-cut on the Cepheid sample and the 
choice of the extincti on curve. The extinction curve of 
IFitzpatrick fc Massal (120071) yields the smallest increase 
in scatter as the colour cut is increased. 


ElV - 7) Max 

CCM89 

F99 

FM07 

0.39 

0.1687 

0.1676 

0.1678 

0.99 

0.1844 

0.1839 

0.1766 


track the scatter in the relations rather than the shift in 
distance modulus, for which the re sults are listed in Ta ble 
2] Assuming the extinction curve of lCardelli et al.l (119891 ). I 
find that the scatter increases from a = 0.1687 at E(V — 
7)Max = 0.39 to a = 0.1844 at E(V — I ) M ax = 0.99. With the 
extinction curve of lFitzpatrickl dl999l ). the scatter increases 
from a = 0.1676 to a = 0.1839. Finally, with the extinction 
curve of lFitzpatrick fc Massal (|2007l ). the scatter increases 
from a = 0.1678 to u = 0.1766. The increase in variance 
for the shallowest extinction curve as the reddening cut is 
increased is only ~ 10 %, less than half what it is for the other 
two extinction curves. 


4 DISCUSSION AND CONCLUSION 

In this investigation I have demonstrated that the Cepheid- 
distance to the spiral galaxy M101 is measured with 
a reduced scatter and reduced sensitivity to systemat- 
ics if I assume the mean Galactic extinction curve of 


































































































6 Nataf 


iFitzpatrick fe Massal 120071) . which is stee per than the 
canon ical Ry = 3.1 extinction curve from ICardelli et al.l 
( 1989 ). The steeper extinction curve increases the distance 
modulus to M101 by ~0.06 mag. The combination of the 
different extinction curve, and the downward adjustment 
in the assumed metallicity of the LMC, yields a distance 
modulus estimate t hat is ~0.09 mag greater than that of 
IShappee fe Stanek d201 ll) . to A/tlmc ~ 10.72 ± 0.03. The 
best-fit metallicity-dependence is dMi/d[ O/H] ss (—0.38 ± 
0.14) mag dex -1 . 

The mean Galac tic extinction curve of 
iFitzpatrick fe Massal (120071 ) yields a superior fit to the 
M101 Cepheids, but it should not be expected to apply 
in general. dKriek fe Conrxnl 120131 ') have shown that the 
attenuation curves of galaxies correlates with their spectral 
type and specific-star-formation rate. Meanwhile, it’s by 
now well-documented that the shape of extinction curve 
is a function o f locat io n even in on e part ic ular galaxy 
dZasowski et al.l 120091 : iPhillips et aid 1 20131 : iDong et al] 
l2014h . I thus argue that a new phenomenological model to 
deal with interstellar extinction is needed, both to under¬ 
stand extragalactic stellar populations and to precisely and 
accurately calibrate the distance ladder. 

The referee brought up the possibility of varying Ry. 
Given that the investigation of this work had access to data 
in only two bandpasses, varying Ry from a fixed family of ex¬ 
tinction curves is indistinguishable from varyin g th e extinc¬ 
tion c urve at fixed Ry. From the formalism of lCardelli et all 
1 19891) . one obtains Ai/E(V — I) = 1.145 at Ry as 2.33 , 
a value common in Type la SNe cosmology ((Phillips et al.l 
120131) . The distinction cannot be visited in this work, it could 
in principle be investigated if there were a larger number of 
bandpasses available, in particular B-band would obviously 
be a significant degeneracy-breaker. 

Some readers may worry that the era of worrying about 
extinction as a so urce o f error in_the distance ladder is end¬ 
ing. For example, iRiess et al.1 (j2012j ) calibrated the near-IR 
period-luminosity relations for a sample of 68 Cepheids in 
the Andromeda galaxy (M31). This strategy is promising, 
as the extinction in WFC3 filter F160W is only ~30% that 
of the extinction in the Landolt /-band filter, and the dis¬ 
crepancy in the predicted value of Afwow between different 
extinction curves is only ~5% of the size of A 5500 ^ (see Ta- 
ble[lj. However, it will also be harder to measure the Hubble 
flow to longer distances given the ~300% reduction in reso¬ 
lution at the longer wavelength, in addition to the possibility 
that Cepheids may be further blended at these longer wave¬ 
lengths, for example by AGB stars which are intrinsically 
bright, red, and numerous. Further, regardless of this issue, 
the existing trove of optical Cepheid data to nearby galax¬ 
ies yields a probe of the interstellar medium that is distinct 
from others currently in use. 


ACKNOWLEDGMENTS 

I thank Benjamin J. Shappee, Brian P. Schmidt, and Mar¬ 
tin Asplund, for helpful discussions. I thank the anonymous 
referee for a constructive review of the manuscript. 

DMN was supported by the Australian Research Council 
grant FL110100012. 


REFERENCES 

Berdnikov L. N., Vozyakova O. V., Dambis A. K., 1996, 
Astronomy Letters, 22, 334 
Bresolin F., 2011, ApJ, 729, 56 

Cardclli J. A., Clayton G. C., Mathis J. S., 1989, ApJ, 345, 
245 

Dong H. et al., 2014, ApJ, 785, 136 

Efstathiou G., 2014, MNRAS, 440, 1138 

Fausnaugh M. M., Kochanek C. S., Gerke J. R., Macri 

L. M., Riess A. G., Stanek K. Z., 2014, ArXiv e-prints 
Fitzpatrick E. L., 1999, PASP, 111, 63 

Fitzpatrick E. L., Massa D., 2007, ApJ, 663, 320 
Freedman W. L. et al., 2001, ApJ, 553, 47 
Gerke J. R., Kochanek C. S., Prieto J. L., Stanek K. Z., 
Macri L. M., 2011, ApJ, 743, 176 
Kriek M., Conroy C., 2013, ApJ, 775, L16 
Landolt A. U., 1992, AJ, 104, 340 

Macri L. M., Stanek K. Z., Bersier D., Greenhill L. J., Reid 

M. J., 2006, ApJ, 652, 1133 
Madore B. F., 1982, ApJ, 253, 575 

Majaess D. J., Turner D. G., Lane D. J., Henden A. A., 
Krajci T., 2011, Journal of the American Association of 
Variable Star Observers (JAAVSO), 39, 122 
Nataf D. M. et al., 2013, ApJ, 769, 88 
Ngeow C. C., 2012, ApJ, 747, 50 

Nishiyama S., Tamura M., Hatano H., Kato D., Tanabe T., 
Sugitani K., Nagata T., 2009, ApJ, 696, 1407 
Paczynski B., Stanek K. Z., Udalski A., Szymanski M., 
Kaluzny J., Kubiak M., Mateo M., 1994, AJ, 107, 2060 
Pejcha O., Kochanek C. S., 2012, ApJ, 748, 107 
Pejcha O., Prieto J. L., 2014, ArXiv e-prints 
Phillips M. M. et al., 2013, ApJ, 779, 38 
Pietrukowicz P. et al., 2012, ApJ, 750, 169 
Planck Collaboration et al., 2013, ArXiv e-prints 
Riess A. G., Fliri J., Valls-Gabaud D., 2012, ApJ, 745, 156 
Riess A. G. et al., 2011, ApJ, 730, 119 
Russell S. C., Dopita M. A., 1990, ApJS, 74, 93 
Sakai S., Ferrarese L., Kennicutt Jr. R. C., Saha A., 2004, 
ApJ, 608, 42 

Schlafly E. F., Finkbeiner D. P., 2011, ApJ, 737, 103 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Shappee B. J., Stanek K. Z., 2011, ApJ, 733, 124 
Spergel D. N. et al., 2003, ApJS, 148, 175 
Storm J., Gieren W., Fouque P., Barnes T. G., Soszynski 
I., Pietrzynski G., Nardetto N., Queloz D., 2011, A&A, 
534, A95 

Sung H., Bessell M. S., 2014, in H.W. Lee, Y.W. Kang, 
K.C. Leung, eds, Astronomical Society of the Pacific Con¬ 
ference Series. Astronomical Society of the Pacific Confer¬ 
ence Series, Vol. 482, pp. 275-277 
Szymanski M. K., Udalski A., Soszynski I., Kubiak M., 
Pietrzynski G., Poleski R., Wyrzykowski L., Ulaczyk K., 
2011, Acta. Astronom., 61, 83 
Udalski A., 2003, ApJ, 590, 284 

Udalski A., Soszynski I., Szymanski M., Kubiak M., 
Pietrzynski G., Wozniak P., Zebrun K., 1999, Acta. As¬ 
tronom., 49, 223 

Zasowski G. et al., 2009, ApJ, 707, 510 












































